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We study quantum compression and decompression of light pulses that carry quantum information 
using a photon-echo quantum memory technique with controllable inhomogeneous broadening of an 
isolated atomic absorption line. We investigate media with differently broadened absorption profiles, 
transverse and longitudinal, finding that the recall efficiency can be as large as unity and that the 
quantum information encoded into the photonic qubits can remain unperturbed. Our results provide 
new insight into reversible light-atom interaction, and are interesting in view of future quantum 
communication networks, where pulse compression and decompression may play an important role 
to increase the qubit rate, or to map quantum information from photonic carriers with large optical 
bandwidth into atomic memories with smaller bandwidth. 
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I. INTRODUCTION 

As any communication system, quantum communication relies on preparing carriers of (quantum) information, 
transmitting those carriers in a reliable way, and processing the information. Obviously, in order to maximize the 
communication rate, the duration of the carriers, naturally photons, should be as short as possible. However, while 
quantum communication channels such as optical fibres or free-space channels allow transmission of broadband, sub 
nanosecond photons, it is often impossible to generate or process such carriers. Of particular concern are key elements 
for quantum repeaters jT] such as certain entangled photon pair sources or quantum memories [5], whose 

bandwidth is often limited by material constraints. 

In this article, we study temporal/bandwidth conversion as a quantum information preserving transformation for 
photonic information carriers. On the one hand, increasing the optical bandwidth in an efficient way, i.e. decreasing 
the duration of an information carrier, will allow increasing the transmission rate over a single quantum channel when 
time-multiplexing several small bandwidth photon sources or photon pair sources. On the other hand, decreasing 
the bandwidth will allow mapping of large-bandwidth photonic carriers into quantum memories with limited spectral 
width. In view of this transformation, we investigate a photon-echo type quantum memory approach based on 
controlled reversible inhomogeneous broadening (CRIB) of a large ensemble of atomic absorbers [SHlOj. Relaxing 
the requirement of symmetric inversion of atomic dctunings A — > — A during absorption of the light and recall, 
respectively, in a more general version with non-symmetric inversion A — > —77 A with compression factor 77 ^ 1, 
we find accelerated/decelerated rephasing of atomic coherences and thus temporal compression/decompression of the 
reemitted light field. We analytically analyze the proposed scheme in optically thick atomic media with transverse and 
longitudinal inhomogeneous broadenings for recall efficiency, fidelity, and gain of the transmission rate over a single 
channel through multiplexing. In transverse inhomogeneously broadened media, the atomic resonance frequencies 
vary normal (transverse) to the spatial coordinate z measured along the propagation direction of the light, and the 
absorption profile is independent of z. In longitudinal inhomogeneously broadened media, the resonance frequencies 
depend linearly on z. Assuming large optical depth for storage and retrieval, we find, for the case of transverse 
broadening, that the recall efficiency is limited by the compression factor 77, while it reaches unity in the case of 
a longitudinal broadened medium. We also find, for transverse broadening, that the fidelity of a recalled photonic 
time-bin qubit with the original qubit is one, regardless the compression factor, but that it is limited in the case of 
longitudinal broadening. We point out that optical pulse compression has previously been considered using traditional 
photon echos and chirped excitation pulses [TTlJ13| . However, similar to data storage [21 [15], this approach is not 
suitable for temporal compression of quantum data. Quantum compression using CRIB has first been discussed in 
[16j . and first observations as well as numerical studies for the case of longitudinal broadening have recently been 
reported [T?|ITg]. 

This article is organized in the following way. We will first discuss the figures of merit chosen to assess the 
performance of quantum compression. We will then introduce photonic time-bin qubits and present two ways to 
describe this basic unit of quantum information. This part is followed by a description of the standard CRIB-based 
quantum memory protocol, which employs a hidden symmetry in the equations describing the atom-light interaction 
during storage and recall. In order to assess the change in the recalled photonic quantum state for CRIB-based 
quantum compression, which does not involve time-reversed quantum dynamics, we have to solve the equations of 
motion. This is done in the next section, where we also derive the efficiency, fidelity and gain of quantum compression 
for the examples of transverse and longitudinal broadening. This is the main part of this article. It is followed by a 
brief feasibility study of our protocol in rare-earth- ion doped crystals. The article terminates with a conclusion. 

II. FIGURES OF MERIT 

Temporal compression of quantum data is of importance to quantum communication, similar to temporal com- 
pression of classical data and classical (tele) communication. However, the criteria imposed on a "good" compression 
procedure are much more severe in the quantum case: reduced efficiency impacts irreversibly on the quantum infor- 
mation rate through photon loss, in opposition to amplitude loss in the classical case, which can be compensated 
by means of optical amplifiers (note that amplification is unsuitable for quantum communication, as stated in the 
no-cloning theorem [19, 20 ). Furthermore, unpredictable modification of the input photonic quantum state \ tpi n ) p 
during compression results in an increased quantum bit error rate (QBER), or requires compensation via not-yet- 
practical quantum error correction |21[ I22j. In opposition, classical information, due to its digital nature, is much 
more tolerant to noise. In this article we therefore use the efficiency e, the fidelity F (which specifies the unpredictable 
change of an input quantum state), and the gain G (which derives from the efficiency and the compression parameter), 
as figures of merit to analyze the performance of quantum compression. 

In the following, we assume pure photonic qubit input states. The compression operation yields a (possibly mixed) 
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photonic output state characterized by a generally not normalized density matrix p ou t,p- We define the efficiency e of 
the compression as 

e = tr(p out , p ) (1) 

The index "p" denotes "photon" and will be added henceforth to avoid confusion of photonic with atomic states. 
Furthermore, we use the following definition of the fidelity F : 

F= p {^in\p'out,p\^in)p (2) 

where p' out p is the renormalized and unitarily transformed density matrix 

P'out.,p=\ U Pout,pU + , (3) 

and p ut,p describes the recalled and compressed photonic qubit. It is obtained from the total density matrix by 
tracing over the degrees of freedom related to the atomic system and possibly non-compressed photonic modes. This 
procedure is justified as one can experimentally restrict photon detection to the desired cases. 

Note that we allowed for deterministic unitary operations U composed of rotations around a x , <r v , or o~ z in the 
qubit Hilbert space to maximize the fidelity. This is similar to quantum teleportation [23], where, depending on the 
result of the Bell state measurement, bit flip, phase flip, or bit and phase flip operations have to be applied to the 
teleported state to recover the initial state. 

Finally, in order to characterize the usefulness of quantum compression in view of enhanced qubit transmission rate, 
we define the gain as 

G = ery. (4) 

This reflects that the detection rate of qubits increases linearly with both the compression efficiency e as well as the 
factor rj by which a qubit can be compressed, i.e. the number of qubits per unit time. 



III. PHOTONIC TIME-BIN QUBITS 



We are concerned with encoding of quantum information into qubits, i.e. quantum states that are described by 

|V)=a|0}+e^|l). (5) 

The coefficients a, (3 and cj> are real, a 2 +/3 2 = 1, and the kets | ) and 1 1 ) form an orthogonal basis in a two-dimensional 
Hilbert space ((i \ j) — 5ij,i,j = [0, 1]). 

In this article, we are specifically interested in so-called time-bin qubits, where the basis states | 0) and 1 1 ) describe 
photon wavcpackcts localized at a particular position z at early and late times t — z/c and t' = z/c + t Qi respectively. 
Differently stated, at a given time t, the photon is in a superposition of being at positions z and z' = z — z a with 
zo = CT o, see Fig. [T] We will use these two pictures interchangeably in our analysis. Time-bin qubits have been shown 
to be very well suited for quantum communication over telecommunication fibres 24, 25]. In the following, we will 
derive a description of time-bin qubits using a physical representation of the abstract qubit Hilbert space. We will 
limit ourselves to one polarization mode. 

A general photonic wavepacket is described in momentum space as [261 127] 

/oo 
dkg(k,t)&+\Q), (6) 
-oo 

where (d^) denotes the creation (annihilation) operator for a photon with wave vector k, [a^, a^,] = S(k — k'), and 
| ) is the vacuum state of light. For a more suitable description of our qubit states, we define operators creating 
forward and backward propagating photons at a particular position (we will index these two cases by / and +, or b 
and — , respectively): 



a+iz) = di (z) = \/l/2ir / dk exp{— ikz}d^, (7) 

Jo 

d~t(z) = d^(z) = \J\j1-K I dk exp{—ikz}d^ . (8) 

J — oo 




FIG. 1. Schematical representation of a photonic time-bin qubit \ tp) = a\0) + f3e"^\ 1), where a and /3 are the amplitudes 
of photon wave packets propagating along the +z axis with relative phase (f>; Sz — cSt, and z = ct . St and r Q denote the 
temporal duration of the wave packets and the relative time delay, respectively. 



The associated annihilation operators are denned by df(z) — (d^(z)) + and d b (z) = (d^(z)) + , satisfying the usual 
bosonic commutation relations. With these new operators, the basis states of a time-bin qubit can be expressed via 
the center z — act and z — z D = ac{t — T ) (a = [+, — ]) of the wavepackets, respectively, their extension 5z and carrier 
frequency ui a : 



(9) 



0) CT ,p = / dz' exp{— iuj a {t — crz'/c)}g(ct — az', 5z)a£(z')\ ), 

J — oo 

= 1 1, Sz, u a ) CTiP 

dz' exp{— iuj a (t — t — az' /c)}g(ct — cr Q — az' , 6z)a^(z') \ ), 



(10) 



= | t — T , 5z,Wo- }o-,p 

where g{ct — z' ', Sz) is the normalized envelope of the photon wave packet with length 8z, which is related to its spectral 
function g(k, 5k) by Fourier transformation: <?(z, Sz) = y/lJ^/K dkg{k, 8k) exp{— ikz}. Note that the identification 
of the states 1 1, Sz, cu a ) a p and 1 1 — t q , Sz, cu a ) ap with | ) a and 1 1 respectively, requires Sz z a = ct so that the 
states are orthogonal and normalized. Using this notation, the general photonic wavepacket can then be expressed as 



\^(t))p = \^(t))f,p+\m)b, P 

= v/l/ 27r X! / d ^<T( t ^')exp{-ia; CT (i- ( Tz7c)}a+(z')|0). 
In the case of the initial, forward propagating time-bin qubit, the amplitude A + takes on the form 



(11) 



A + {t —} — oo, z) = V27r{ug(ct — z, Sz) + exp{iw + r + i(/)}(3g(ct — z — ct , Sz)} 
and A_{t — > — oo, z) — 0, hence 



(12) 



I ipin(t) )p = a\ t, Sz, lj + ) fiP + e^Pl t - t , Sz, uj + ) ftP . 



(13) 



A schematical representation of a time-bin qubit assuming Gaussian shapes of the basis wave packets is given in Fig. 

m 



IV. QUANTUM STATE STORAGE BASED ON SYMMETRIC CRIB 

An extensively investigated approach to quantum state storage is based on controlled reversible inhomogeneous 
broadening (CRIB) of an isolated absorption line [BHTU], It has recently led to storage efficiencies up to 69% [25], 
and the possibility for random access quantum memory has been demonstrated [17] . In the ideal, standard CRIB 
scheme, the optical input data is launched along the forward (+z) direction into a resonant, optically thick atomic 
medium (olL —¥ oo, where a is the resonant absorption coefficient and L is the length of the atomic medium). The 
spectral width Sui n of the input light should be narrow in comparison with the inhomogeneous broadening Aj n /j of 
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the resonant atomic transition lf)2, i.e. Suji n < Aj n /j, and larger than the homogeneous line width j eg of the optical 
transition, i.e. 6u>i n > ^ eg . 

All atoms are initially prepared in the pure state \g) = Y\.j=\®\9)j (where \g)j is a long-lived (ground) state of 
the j-th atom). Thus the initial light-atom state, denoted using the capital letter (for t — > — oo) is 



| *(t ->■ -oo) ) f = | ^, n (t) )/, p (8) | .9 ) . (14) 

After absorption of the light, the atomic coherences acquire phase factors exp{— i(Aj + w eg )i}, where is the 
detuning of atom j as compared to the central absorption frequency uj egi leading to fast dephasing of the excited 
collective atomic coherence. The absorption process can be time reversed by inverting the atomic detunings at time 
t\ (i.e. Aj — > —Aj) and applying a position dependent phase shift of 2kz. This leads to rephasing of the collective 
atomic coherence, hence irradiation of an echo signal at time t ec ho in backwards direction as a perfect time reversed 
copy of the input data jB] . This reversibility is based on a hidden symmetry of the equations describing the evolution 
of the slowly varying light and atomic parameters during storage and retrieval [9]. 

Analyzing the evolution of the complete wavefunction (including the fast, time varying part exp{—ick eg t} where ck eg 
is the central frequency of the atomic transition), it was shown [29) that the photon operators (in momentum space) 
of the input (forward) and output (backward) light fields are related by exp{— ic(k + k eg )t}a k+k — > — exp{ic(k — 

)}al_ k . As we will show in more detail in Appendix A, this results for input photonic time-bin qubit 



states (Eq. 13) to emerge as 



I 4>out(t - t echo ) ) btP = - (a\ t, 8z,u- ) biP + e ^+ 2 <w°) f3\ t + r , Sz,uj^ ) 6)P J , 



(15) 



where we have ignored a global phase shift, and where w_ = 2uj eg — uj + denotes the new carrier frequency. Hence, 
the recalled qubit state is associated with an exchange of the leading and trailing wavepackets, and the now leading 
wavepacket (with probability amplitude proportional to j3) has acquired an additional phase that depends on the 
time delay r Q between the wavepackets, and the atomic transition frequency u> eg (see also [30 ). Please note that 
not only the order of the wavepackets change, but that each wavepacket's temporal envelope is also time-reversed. 
For simplicity of notation, we will henceforth restrict our investigation to time-symmetric, Gaussian shaped basis 
wavepackets described by g(ct — z, 6z) = ^/2/ (tt(Sz) 2 ) exp{— (ct/Sz) 2 }, where t = t — z/c. 

Returning to the abstract qubit notation, we find that the initial state | ipi n ) = a\ 0) + e**/3| 1 ) is transformed in 
the quantum memory into | ipout ) = cx| 1) + e l ^ +2u) ° aT °^ /3| 0), i.e. that the two states are related by a deterministic 
unitary transformation T 

IVw) =T\^ m ) =e tu ^°a x R z (6)\^n), (16) 

.0 

where R z (9) = e" 1 z denotes a rotation of 9 — 2uj eg T around a z , and a x is the bit flip operator. 

Hence, for the ideal standard CRIB protocol, we find the efficiency to be one. Furthermore, from Eq. we find 
U = T~ l , hence p' out p = Pi n , P and thus F=l. Please note that two subsequent storage sequences lead to compensation 
of the bit flip as well as the additional phase factor of 2ui eg T . 

Finally, being obvious in the case of "no compression" and unity efficiency, we note that the gain G is one. 



V. QUANTUM COMPRESSION BASED ON GENERALIZED CRIB 



In the CRIB protocol described above, the efficiency, fidelity, and gain can be derived using arguments stemming 
from symmetries in the equations of motion. Relaxing the perfect reversibility of atomic detunings by introducing a 
more general relation between the initial (t <t\) and output (t > t\) spectral detunings 

A j (t>t 1 ) = -TjA j (t<t 1 ), (17) 

with rj being the compression parameter, this symmetry-based approach is no longer possible. Here, we pursue the 
following approach: We start with a general photonic qubit state (Eq. [5| described in the abstract two-dimensional 
Hilbert space. We then express the qubit state using a physical description in position space (Eqs. ([9 10), and (13), 



respectively). Using the equations of motion, which we will introduce in the following section, we then calculate the 



state of the combined atom-photon system after quantum compression. From Eq. (15), we find the single photon 
density matrix. Redefining new basis states | ) and | 1 ) , determined by the compression parameter, we then express 
the output state in qubit state notation, which finally allows us to calculate our figures of merit (Eqs. (JlJ, pj), and 
Q, and the unitary operation U (see Eq. ([3])). 
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A. Basic equations 

The interaction of the quantum fields with atoms is given by the Hamiltonian 26, 271 [29] 



H(t)=H +Hx{t), (18) 

where 

JV 

(19) 



and 



Hi(t) = 

N 



hgV2nY^ \ {%(%) + a>b(zj)}PZ g + h, 

N 



3=1 

/OO Q 
dza+(z) — a a (z). (20) 
-oo OZ 

Hq describes the total number of excitations in the system and commutes with the total Hamiltonian H. The first 
term in H\ (t) describes the atom-field interaction with atomic operator P^ n = | m ) jj { n \ . N denotes the number of 
atoms, g{oS) = g = id ge ( 2 ^ g g ) 1 ^ 2 , di2 is the dipole moment of the atomic transition \g) O |e) , S is the cross section of 
the light beams, e is the electric permittivity, and h is Planck's constant divided by 2ir. The second term in Eq. [20] 
describes inhomogeneous and homogeneous atomic line broadening within a unitary approach to quantum evolution. 
Indeed, it has been shown [311 E2] that the decay rate of the atomic coherence 7 eg can be calculated by a statistical 
average over local stochastic phase fluctuations S<fii g (t,t') = f., dt n 8E^ g (t)/h of the atomic transition: 

( exp{-z^ eg (t, f )}) - cxp{- 7e3 (i - t')}, (21) 

where (6^ eg {t, t')} — 0. These fluctuations are due to interaction between the absorbers and its environment. Finally, 

the third term of Hi (t) takes into account the spectrum of the localized photonic wave packets (here expressed through 
position operators.) 



Using the total Hamiltonian (18), the initial photon-atom state (14) will evolve into 



!*(*)) = l*(*))p + l*(*))a» ( 22 ) 



l*(*)> P = M*)>P®lff>, (23) 

(with | ip(t) ) p as defined in Eq. (|6])) and 

N 

!*(*)>„ = E 6 ^4IS>® I 0>> (24) 

3=1 

with bj(t) being the probability amplitude for atom j to be in the excited state. Initially bj(t — > — oo) = 0, i.e. all 
atoms arc in the ground state, and A-(t, — > — oo, z)) = 0, i.e. all backward modes are empty. 



In the Schrodinger equation with the Hamiltonian introduced in (Eq. 18), the quantum evolution during absorption 
of the forward propagating light field (a = + for t < ti) and for re-emission of the backward propagating light field 
[a = — for t > ti) is given by 
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TIME: 

a) t<0 

b) t = 

c) t = ti 

d) t=(1+1/r|)ti 

e) t>(1+1/r|)ti 



P a 

it 



Atomic medium 



5z 




t 



FIG. 2. Schematics showing different instances in the quantum compression protocol in the case of transverse broadening. The 
figure depicts forwards and backwards propagating photonic wave packets, and atomic excitation. 



'1JL 

v c dt 



dz 



) A a (t, z) = i(Trn Sg* /c) exp{iw cr (t — az/c)}b(t, z) 



si 

dt u 3 



bj(t) = -i{uj eg + &j(t) + SAi (t)}bj(t) + igA a {t, Zj ) exp{-iw CT (t - aZj/c)}, 



(25) 
(26) 



where b(t, z) = (n c 5) _1 X^=i bj(t)5(z — zj) is a collective atomic variable describing the averaged, position dependent 
atomic coherence, and n is the atomic density. Using these equations, we will now evaluate quantum compression of 
photonic time-bin qubits in atomic systems with transverse and longitudinal inhomogeneous broadening. 



B. Transverse broadened media 



Transverse broadening is the inhomogeneous broadening occuring naturally. In such media, the atomic absorption 
profile is independent of the spatial coordinate z measured along the propagation direction of the light. In rare-earth- 
ion doped inorganic crystals (RE crystals) with low site symmetry and sufficiently narrow absorption line, controlled 
transverse broadening can be introduced through the Stark effect [33] by applying an electric field gradient transverse 
to the propagation direction of light. For RE ions implemented in a glassy host, e.g. an optical fibre, controlled 
transverse broadening is already obtained for uniform electric fields, due to the random orientation and magnitude of 
the electric dipole moments of the RE ions [31] . In the following, we assume a broadened medium with Lorentzian 
lincshape G(A/Aj„^) = Ai„f l /((Af nh + A 2 )tt) with inhomogeneous spectral width A„^ 3> cSk. The center of the first 
wave packet |0)/ jP enters the atomic medium at t = (see Fig. ph. During the interaction with the medium, the 
light field is partially or completely absorbed, depending on the optical depth. At time t\, after sufficient dephasing, 
we change the detuning A 3 to — ^Aj, and we also apply a position dependent phase shift of (2k eg — Sk)z, which allows 
phase matching for the retrieval of the light field in backwards direction [8, 9, 35, 36J. 5k describes a small deviation 
from the perfect phase matching. This leads to reemission of the light field, with the center of the (now trailing) 
wavepacket exiting the medium at time t — t ec } lo . The atom-light state for t 3> t ec h Q is then given by 



I Ht » t echo ) } = | Ipfit) ) p ® | ) + | *(t) )a + I M*) >p ® I 0>- (27) 

As visualized in Figs. [5J:, d and e, the first term describes the damped, non-absorbed photonic component \ipf(t))p = 
exp{—a L/2}\ ?pi n (t) ) p that continued to travel in forward direction behind the atomic medium, the second term the 
remaining atomic excitation, and the third term denotes the now backwards moving, retrieved photon field, which is 
the subject of our investigation. a is the on-resonant absorption coefficient, L the length of the medium, and a a L 



the optical depth. As shown in Appendix |VII A the amplitude of the retrieved photon field, assuming for simplicity 



a Gaussian spectral shape, is described in abstract qubit notation by: 



| Mt) ><*> = R (t) \ e^ 1+1 ^^a\ t, Sz',u>- ) b + e**'p\ t + r' , 6z' , w_ > 6 1, (28) 



where the superscript "(t)" denotes transverse broadening, where we ignored a global phase shift, and where <p' 
4> + (1 + l/rj)uj e gT . Furthermore, 



/CO 
dz' exp{-ioj a {t^ - t echo + z'/c)}g{ct^ - ct echo + z' , 5z')a+(z')\ ), 
-OO 



(29) 



where denotes the early or late wave packet (t^ m ' € [t,t + T ]), techo = (l + l/»7)ti — 6tu, 6tn = i^ eg [l + r]](5z' /c) 2 , 
t' d = T /rj and 6z' = 6zjr\. The factor 

= V leg {t 1 )M^{5k,a L)e 1 J 2 {r 1 ) (30) 

combines different sources that affect the efficiency of the quantum compression through phase relaxation 7 eg , phase 
mismatch 5k, and limited optical depth a a L. Its upper limit, given by e D (not to be confused with the electric 
permittivity e ), is determined by the compression parameter r\. Specifically, we find: 

IV (i x ) = exp{-(l + l/vhe g {h -t - \r]6t R )}, (31) 



M(*) (6k, a L) = l ~ e ^t^ lh) ?r L \l 5kL} - (32) 



Let us now relate Eq. ( [28] ) with our figures of merit. First, it is important to note that Eq. (28) describes again 
a photonic time-bin qubit with modified separation and width of basis states. Indeed, comparing the new with the 
old basis states (Eq. ([13])), we find that the separation between the two basis wavepackets has changed from z a 
to z' a — z /rj, and their extension from Sz to Sz' = Sz/rj, while the wave packet amplitude has changed oc ry 1 / 2 as 
visualized in Figs. [2]d and e. 

Second, it is important to realize that the quantum information encoded into a qubit is independent of the physical 
realization of its abstract basis states |0) and 1 1). This allows us to relabel the basis states after compression, i.e. 
1 1, Sz', w_ ) — > | 1 ), and 1 1 + t' d , Sz' . w_ ) — > | ), resulting in: 



(34) 



Now, taking the trace of p ou t, P = \ V 1 )pp{ i> I we nn d the efficiency according to Eq. |I]) 

e W = (r 7c3 (t 1 )) 2 |M(*)(^,a L)| 2 e (r,){a 2 e- 2 ( 1+1 /"^^ +/3 2 }. (35) 

We now discuss our figures of merit for some particular cases. Ignoring phase relaxation (i.e. 7 e9 = 0) and assuming 
perfect phase matching (i.e. 6k = 0) we find the efficiency to be given by: 

eW = (^TIF (1 " exp{ ^(! + iA?HL}) 2 . (36) 

It is depicted in Figure|3]as a function of the optical depth a^L and the compression parameter r\. First, we note that 

for standard CRIB (without compression, i.e. r\ — 1), we find the previously published result = (I — exp{— a L}) 2 
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FIG. 3. Recall efficiency for transverse broadening as a function of the compression parameter r\ and the optical depth a a L 
(see Eq. [Esg}). 



[551 [37]. Second, we point out that for a given compression factor, the maximum efficiency is always obtained for 
infinite optical depth, and hence the upper limit e max is only determined by the compression factor: 



At) 



(V) 



4;; 



(1 + riY 



(37) 



To give an example, a quantum compression with 3 — 2\/2 < i] < 3 + 2\/2 results in a recall efficiency limited to 50%. 

To assess the fidelity, we compare Eq. (34) with the initial qubit state in Eq. We find that the re-normalized, 
compressed qubit state | tf> ) p can be unitarily transformed into the initial state by a o~ z rotation of angle $ — cf> = 
(1 + l/r])ui e gT and a bit flip operation o x . Furthermore, one should take into account an additional o~ x rotation to 
compensate for the amplitude reduction factor exp{— (1 + 1 /r])"f e gT } arising in the case of significant phase relaxation 
during time (1 + \ jr\)T . This immediately implies that the fidelity as defined in Eq. ^ is always one, regardless the 
compression factor, phase relaxation, or phase mismatch. Note that the fact that atomic phase relaxation does not 
affect the fidelity in time qubit storage has been inferred from stimulated photon echo experiments with intense light 
pulses [38 . 

Figure [4] depicts the gain as a function of optical depth and compression parameter for transverse broadened media. 
As an example, assuming a a L=2, we find an increased communication rate, i.e. >1, for r\ > 1.7. For large 

(infinite) optical depth (i.e. maximum efficiency as described by Eq. (37)), the gain is upper bounded by 



G {t) {aL -> 00) = 4?7 2 /(l + rff\ 



»7»1 



(38) 



a modest, yet significant improvement over quantum communication schemes without compression. 

To finish this section, let us briefly inspect Eqs. (30)-( [33| in view of symmetry between compression and decom- 
pression, i.e. under exchange of rj to I/77. First, we recall that the fidelity is one, regardless the compression factor: 
F(j]) = F{l/rf) = 1. Second, we find that the upper limit of the efficiency is symmetric with respect to compres- 
sion/decompression: Cmaxiji) = e max{^-M) < 1- However, limited optical depth, phase mismatch, and atomic phase 
relaxation break the symmetry in the efficiency. It should be possible to demonstrate this surprising behavior using 
only symmetry arguments and the equations of motion describing the recall efficiency, without having to solve the 
equations. This is likely to lead to a more profound understanding of the physical principles of the here studied 
compression/decompression protocol. 



C. Longitudinally broadened media 



Longitudinal broadening is an inhomogeneous broadening that can not be found naturally. In this case, for each 
position z in the medium, the atomic absorption profile is given by a narrow line (here for simplicity assumed to be 
homogeneously broadenend), whose center frequency varies monotonously with z. Longitudinal broadening can be 
realized through the Stark effect in non centro-symmetric RE ion doped inorganic crystals by applying a electric field 
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FIG. 4. Gain for transverse broadening as a function of the compression parameter r\ and the optical depth a L. 

gradient longitudinal to the propagation direction of light [7|. In the following, we assume absorption lines whose 
detuning with respect to the light carrier frequency varies linearly with position z within the crystal, which extends 
from z = -L/2 to z = +L/2: 



-XZ- 



(39) 



The initial inhomogeneous absorption line, obtained after integration over all atomic positions z, is assumed to be 
broad compared to the spectrum of the photonic wavepacket: Ai n h — \L 3> c8k. A schematics of the compression 
procedure depicting relevant instances of the protocol is given in Fig. [5j 

A similar approach as in the case of transverse broadening leads to the following wavefunction for the backwards 
emitted light field (for more details see Appendix VII B ) : 



| Mt) ><° = R (l) \ e-( 1+1 ^h^° a \ t, t - 4k. h + t, t - tf k + t' q , 8z', ^ ) b 



(40) 

where the subscript "(1)" denotes longitudinal broadening and where we ignored again a global phase shift. Further- 
more, 

\t,& l \8z',uj_) b = 

/OO 
dz' exp{i${t + z7c,77)}exp{-io;_(t (m) - t echo + z' /c)}g{ct {m) - ct echo + z' ',6z')a£ \z')\$) , 
-OO 

and e[t- tf k +T' ,t- 42] with 



(41) 



4 ) -^sb:(^ + 2(C/x)(i-i/»7)). 



(42) 



The parameters 8z' = 8z/rj, t' — t /t], and <p' are the same as in the transverse case, and £ = im Sg 2 / c. Furthermore, 

^ ) =r 7cs (ti+41 ) )Af (i) (C/x^) (43) 



with 



MW(C/ X ,»7) = v /(l-exp{-^})(l-exp{-^}). 
Finally, the time dependent phase $(r = t + z/c, rj) is given by 



(44) 
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Time: 

a) t<0 

b) t = 

c) t = ti 

d) t=(1+1/Tl)tl+tsk 



(i) 



e) t>(1+1/r|)ti+t 



6k 
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FIG. 5. Schematic showing different instances in the quantum compression protocol in the case of longitudinal broadening. 
The figure depicts forwards and backwards propagating photonic wave packets, and atomic excitation. Note the difference of 
the re-emission time t ec ho and the different localization of atomic excitation compared to the transverse case (see Fig. [2). The 
change of colour of the irradiated wave packets denotes additional frequency shifts (see Eq.(46l). 



where r' m = r m /r) 2 , and r TO = 2 



For time-bin qubits with sufficiently narrow basis wavepackets, this nonlinear 



phase change translates into a different frequency change of each wavepacket (8u)q, 8u>i) as compared to the carrier 
frequency. From 8ui(t) — d<f>(t)/dt we find: 



Suj 



5uji = — 
X 



1 



V 



h-r + r m + 8k/x h-T + T„ 

1 r\ 
h + r m + 8k jx h + t, 



(46) 



Note that a given dephasing time t determines close-to-homogeneously broadened slices of length 81 in the atomic 
medium with 81 « (\t) _1 where coherence, hence radiation, remains. We believe that coupling between light and 
collective atomic coherence (as included in the analyzed equations) within these slices plays an important role for 
those frequency shifts, which are proportional to the number of atoms n — NSl in each slice, and the square of the 
photon-atom coupling constant g (i.e. Sou oc g 2 NSl). The nonlinear phase modulation also leads to an additional, 

possibly substantial phase difference between the two basic wavepackets after recall: Sipoi = f ^ tcc ^° +To /' , dt8uj(t) — 
^{techo + Ta/v) ~ & (techo)- The phase difference is depicted in Figs. [6] and [7j The effects of the non-linear phase 
shift are enhanced with increasing effective optical depth K e ff = 27r£/x, time delay r between the wave packets, and 
diminishes as the dephasing time t\ increases, and the compression parameter 77 approaches unity. 

Phase and frequencies shifts in the here analyzed case of backward emission are a result of the presence of phase 
mismatch (Sk^ ^ 0), or compression/decompression (77 ^ 1), or both, but they do not occur in the phase-matched, 
reversible case. However, in the case of forward emission, the phase and frequency shifts also arise for non-compressed 
recall (77 = 1), provided the rephasing and dephasing times t are bounded [39, 40J. 

Let us now discuss our figures of merit in the case of longitudinal broadening and compare it with transverse 
broadening. As mentioned above, the new basis states (Eqs. (41 1) do not transfer into each other through time 
translation alone, in contrast to the case of transverse broadening. For instance, assuming 8t=100 ns, r D = 200 ns, 
2ttC/x = 6, 77 = 3 and a sufficiently large dephasing time t\/8t = 20, we find a frequency difference | 8ujq — 8u>i | s» 2nx 
17.32 kHz . As this value is small compared to the spectral width of the wave packet c8k s» 2ttx 1.59 MHz, we can 
safely ignore this effect as a limitation to the fidelity. It is thus possible to map the new basis states onto the qubit 
basis states | ) and | 1 ) . Furthermore, as before, we can compensate for the additional relative phase arising during 
compression using an appropriate a z rotation. Hence, the fidelity in longitudinal compression can be close to one. 
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FIG. 6. Phase shift Sipoi for quantum compression as a function of the time delay r and the compression parameter rj, with 
dephasing time ti = 20St, Ai n hSt = 10, and 2n(/x = 67r. St denotes the temporal duration of the initial basis wavepackets, 
and is the effective optical depth. 




FIG. 7. Phase shift Sipoi for quantum decompression as a function of the time delay r D and the compression parameter r\, with 
dephasing time t\ — 205t, A in hSt — 10, and 2-7r£/x = 67r. 



The recall efficiency is given by 



c(0 



(R {1} ) 2 = (T 7 .„(ti +*S)) 3 I Af (0 (C/x,»7) | 2 {a 2 e - 2(1+1/ ^— +/? 2 }- 



(47) 



As compared to Eq. 

storage process (Eq. 
phase mismatch only 



+ ^a2)) 2 j which characterized atomic decay during the 



35), we find the same function (I\ (ti -r u Sk 

31)), only with slightly modified argument (now including an additional delay tgl). Note that 
eads to a decrease of the efficiency when paired with significant atomic phase relaxation 7 eg , in 
contrast to the transverse case (Eq. (|32|). 

Figure [8] depicts the efficiency in the case of negligible phase relaxation (i.e. 7 e9 = 0) as a function of effective 
optical depth 2tt£/x and compression parameter 77. Furthermore, a comparison of efficiencies of the longitudinal and 
transverse compression schemes for various initial optical depths is shown in Fig. [9j For small initial optical depth 
( K eff = 2ir(/x ~ 0.57T, or n e ff = a Q L < 0.57T, respectively) and large compression parameter (K e ff/rj < 1), both 
schemes feature the same efficiency: 



(48) 



This behavior is expected as reabsorption in transverse as well as longitudinally broadened media is negligible in the 
case of small optical depth. However, as the optical depth increases, reabsorption becomes more and more important 
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FIG. 8. Recall efficiency for longitudinal broadening as a function of the compression parameter rj and the effective optical 
depth 2ttC/x (see Eq.(d47b). 



40,(0 




2 4 6 8 10 



FIG. 9. Recall efficiencies for compression schemes employing transverse (red-dashed lines) and longitudinal broadening (blue 
lines). The initial optical depths K e ff = a a L = 2nC,/x are 0.l7r,7r,47r (bottom to top sets of two curves). For s mal l optical 
depth and sufficiently large compression factor rj, the efficiencies e'*'(j?) and e"'(f?) are equal, as described by Eq. (481. As the 
optical depth increases, schemes based on longitudinal broadening perform better. 



in the case of transverse broadening, reflecting the departure from reversible light-atom interaction, while it is no 
issue in the case of longitudinal broadening, due to the correlation between atomic detuning A and position z (see 
Eq. (39)). In particular, for longitudinal broadening, we find 



e (i) (f? ->• 0) = (1 - exp{- Ke// }) = l| Ke// »i, 
while decompression with transverse broadening yields a monotonously decreasing efficiency 



(n -> 0) = 4r ? | Kc//>>1 



(49) 



(50) 



and eventually becomes zero. 

Furthermore, the compression efficiency is only limited by the optical depth re e // during absorption, as well as 
by atomic relaxation. In particular, there is no compression factor dependent upper limit, in opposition to e^' (see 
Eq. (37)). Hence, provided sufficient optical depth 2ttQ/x > 1 an d 2n(/(r]x) > 1 (which may be difficult to achieve 
for large compression parameters), one can always achieve close-to-unit efficiency. This gives longitudinal broadening 
a clear advantage over transverse broadening for quantum compression. 

Obviously, the improved performance is reflected in the gain. As shown in Fig. (10), we find that the gain can 



always be increased when increasing the effective optical depth or the compression parameter, as opposed to the case 
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FIG. 10. Gain for longitudinal broadening as function of compression parameter rj and effective optical depth K e ff = 2tt(/x- 



of transverse broadening where the gain was limited to four (see Eq. (38l). Indeed, using Eqs. Q and (47), and 
ignoring atomic relaxation (i.e. 7 e9 =0), we find 



G (l) = 77(1 - exp{-K ef f)/r)}(l - exp{-K e// }). 
Hence, assuming sufficiently large effective optical depth, we find 



(51) 



G^(r,/K eff ^0)=r 1 \ Keff 



(52) 



VI. FEASIBILITY STUDY 



Before we conclude this article, let us briefly discuss possible implementations of our proposal. Quantum compression 
as considered here requires an atomic ensemble with (at least) three suitable energy levels. Two levels must be coupled 
resonantly by the photonic wavepacket to be compressed. The third level, generally another ground state, is employed 
for temporal mapping of the excited optical coherence onto ground state coherence using two counter-propagating 
7r-pulses. This results in the implementation of the 2kz phase shift, i.e. in emission of light in backward direction [B]. 
The width of the optical absorption line, which may be inhomogencously broadened, should be sufficiently small. It 
should be possible to broaden the line in a controlled way via external fields, and the optical depth after broadening 
should be large. The maximum bandwidth of the optical wavepacket to be absorbed or re-emitted is given by the width 
of the externally broadened atomic transition, while the minimum bandwidth of the optical wavepacket is limited by 
the atomic linewidth before broadening (more precisely, the minumum optical bandwidth is a small multiple of the 
atomic linewidth as some controlled broadening is required for controlled dephasing or rephasing). These conditions 
are similar to those required for quantum state storage based on standard CRIB [TU] , but without the need for a long 
storage time, which considerably relaxes material requirements. 

Any material that fulfills these requirements would be a suitable candidate for quantum compression. We note that 
some aspects of the here presented theoretical study have recently been observed using atomic vapour, longitudinal 
broadening, and bright pulses of light [17j . In the following, we will discuss RE- ion doped solid-state material in view 
of quantum compression. 

RE doped solids have been studied extensively for data storage experiments based on stimulated photon echoes 
[41j , and are currently being investigated for photon-echo quantum memory [5] HOj . When implemented into crystals 
and cooled down to temperatures below 4 Kelvin, RE ions typically feature homogeneous linewidths (for optical 
transitions) of a few kHz, and values as small as 50 Hz have been reported [32]. Linewidths in glasses are larger, but 
can still be around or below 1MHz [43H45] . The transitions in crystals and glasses are inhomogeoeusly broadened, 
with values ranging from 40 MHz [46] in crystals to hundreds of GHz in crystals or glasses j42j [47] . The preparation 
of the initial, narrow absorption line therefore requires an initial preparation step via optical pumping |48H51j . In 
low symmetry hosts, the RE ion's quantum states acquire permanent electric dipole moments, which allows shifting 
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their energy levels through the application of a dc electric field - the dc Stark shift (for a review on optical Stark 
spectroscopy of solids see |33j). A different dipole moment for two energy levels consequently results in a shift of the 
associated transition frequency, allowing, through the application of a field gradient, to broaden an isolated absorption 
line. Typical frequency shifts in RE ion doped materials are ten to hundred kHz for an electric field of one VcmT 1 . 

Praseodymium doped Yttrium Orthosilicate (Pr 3+ Y2SiC>5) is a very promising candidate for the demonstration 
of quantum compression in RE solids. This crystal has been employed for storage of light pulses using longitudinal 
broadening [35], and recall efficiencies up to 69% have recently been reported [28 .Taking into account a spectral width 
of the 3 i?4 <^> 1 £>2 (605.9 nm wavelength) absorption line of 30 kHz after optical pumping [32], and a Stark coefficient 
of 112.1kHz/(Vcm- 1 ) [12], we find that Pr 3+ Y 2 Si0 5 allows absorption and reemission of photonic wavepackets with 
bandwidths between approx. hundred kHz and a few MHz, i.e. temporal durations between approx. hundred nsec 
and a few /xsec. This results in a maximum compression (scaling) factor r\ of ~100 (or rj ~1/100). Taking into account 
the demonstrated high efficiency storage [28], one can expect a gain > 10. However, note that the maximum 
bandwidth for this transition is ~10 MHz, due to limitations imposed by ground state hyperfine splitting on the 
preparation of the initial absorption line via optical pumping [49] . 

For RE materials with larger ground state splitting, allowing storage or recall of shorter wavepackets, it is interesting 
to consider RE-ion doped crystalline and amorphous waveguides. Due to the possibility to implement electrodes with 
spacing as small as ~10 jum, large electric fields, i.e. large Stark shifts, can be obtained through application of modest 
voltages. So far, waveguides in Erbium or Thulium doped Lithiumniobate crystals and Erbium doped silicate fibres 
have been investigated in view of quantum state storage [Ml [35] 2U [53] , and Stark broadening up to ^100 MHz seems 
feasible, allowing for storage or recall of pulses as short as ~10 ns. 

VII. CONCLUSION 

In summary, we have studied quantum compression and decompression of photonic time-bin qubit states employing 
a generalized version of CRIB-based photon-echo quantum memory. Assuming high optical depth for storage and 
retrieval, we find, for the case of transverse broadening, that the recall efficiency is limited by the compression factor 
77, while it reaches unity in the case of a longitudinal broadened medium. We also find, for transverse broadening, that 
the fidelity of recalled photonic time-bin qubits with the original qubit is one, regardless the compression factor, but 
that it is limited in the case of longitudinal broadening. Taking into account realistic experimental data, we foresee 
that quantum compression will be useful for quantum communication and computation applications. In particular, 
it enables enhancing the data rate in quantum communication schemes through temporal multiplexing, and allows 
mapping of broadband photons into small-bandwidth quantum memory. 

Our analysis reveals new aspects of coherent photon-atom interaction, specifically photon-echo type interactions, and 
highlights the advantage of schemes employing longitudinal broadening over transverse broadening. The theory can 
be generalized in a straightforward way to encoding of quantum information into multiple (n > 2) (discrete) temporal 
modes of photons, including compression of intense light fields (provided the number of photons is smaller then 
number of resonant atoms). However, quantum compression/decompression of light carrying quantum information 
with continuous spectrum (continuous quantum variables) requires further theoretical investigations that take into 
account added quantum noise in the case of limited recall efficiency. Finally, we note that pulse compression has been 
observed in the context of frequency conversion based on Raman adiabatic transfer of optical states [54] . a protocol 
based on electromagnetically induced transparency that may provide an interesting alternative to our approach. 
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A. Appendix — Transverse broadening 



First, we find a solution of Eqs. (25) and (26) given an input photon amplitude A+(t < t\,z) and excited atomic 
coherence bj(t < t\). Assuming all atoms to be initially in the ground state (bj(r = t Q — > — oo) — > 0) and using 
a temporal Laplace transformation for the amplitude Ai p . + (z) — j d,TA+(t,z)exjp{—p(t — t a )} we find the formal 
solution of Eq. p6|) 



t oo 

bj(t < ti) = ig— cxp{—iuj + (t — t Q )} J dt' J disA u ^ + (zj)exp{—ii't' — i(u eg + Aj—Lj + )(t — t') — i5(f>i g (t,t')}, (53) 

t Q — OO 

where we have also used a backward Laplace transformation A + (t, z) = (2tt)^ 1 di / A v ^ + (z) exp{— iu(t — t )} with 
p = —iv. Putting Eq. (53) in Eq. (25) and taking into account phase relaxation as described in Eq. |2l"] ) we find the 
following expression for the sum of arbitrary atomic functions Fj(t, zj): 



N 



V F^t, z } ) exp{-i / dt'SAi g (t')}S(z - Zj ) \ N>1 

3=1 Jt ° 



= (n S)exp{- leg (t 1 ~t )} J dAG(A/A mh )F(A,t,z). 

— OO 

The Fourier components of the input field are given by 



(54) 



A v<+ {z) = exp{-ia + (j/)z}A^, + (0), 



where 



(55) 



a + (v) = a (~/ eg + A inh ) 



dA 



G(A/A mh ) 



[leg ~ l{v + 0J + - 0J eg - A)} 



(56) 



is the frequency-dependent absorption coefficient for an arbitrary, inhomogeneously broadened absorption line 
G(A/ A inh ) of the transverse type, and the on-resonant absorption coefficient a a — 2Trn Sg 2 /{c(-f eg + A inh )}. The 
amplitude A + (t > t±, z — L) in Eq. (55) describes the light field at z > L after interaction with the atomic medium. 
The Fourier component at the beginning of the medium ^^(O) is determined by the initial state of the input light 
at z = 0. 

We now consider inhomogeneous broadening with Lorentzian lineshape G(A/Aj„^) = ^a^+a? — )' assumm S the 
spectral width of the input fiel d n arrow compared to the atomic line width (Su> + <S A in /j). This results in a + (v) — 

3 + A,„ h ) 



[7e g +A i „, l +i(a 



in Eq. (56 I 



To allow for light re-emission m backward direction, we change the phases of the atomic amplitudes through the 
application of two short 7r laser pulses that temporally map atomic coherence on the optically excited transition 
onto a ground state transition [6l 135] . as discussed in Sec. |IV| In general, this phase matching operation leads 
to some phase mismatch between the atomic coherence and the irradiated echo signal that we quantify through 



an additional atomic phase shift Skzj: bj(t\ + (e — > 0)) = bj(ti) exp{idkzj} in Eqs. (25) and (26) (with index 
a = — denoting echo signal emission and Aj(t > t\) = —TjAj ). T aking into account the initial atomic amplitudes 
bj{t) — bj(ti) exp{iSkzj + ir]Aj(t — ti) — iS4>l g (t, ti)} and Eq. (54), we find the irradiated light field to be described 
by: 



g z Ai Pt -(z) 



= {7rn Sg 2 /c) / df exp{-pf}P echo (f, z) 



Aip,— (z) 



P + (leg + V^-inh) + i(^eg ~ 



(57) 
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where P ec ho(j, z) is the rephased atomic coherence that serves as a source for the irradiated field: 

Pechoir, z) = exp{-(l + r))j eg (r - h)} 

exp{i[5kz + cj_ - uj eg + rj{u + — uj eg )]f + i(uj eg — w+)(l + 



2Ajnh exp{iv[i](f - ti) - t{\} 
[{uj+ + v-u eg + i leg ) 2 + Af n/t ] L 2[1 + i(u} eg -W+- v)/( leg + A inh )] 



rfl/ 77^ ~ " , JZ. N2 , a 2 i ex P{- n[1 , J7. . " 7X777 vr}^,+ (0). (58) 



We used again the temporal Laplace transformation for light and atom amplitudes within the temporal interval oo], 
new variables f = t + z/c, z = z, and we have taken into account that initially A—(t\, < z < L) = 0. The solution 
of Eq. (57) has the form of a double integral 

A-(f,z) = 

OO _ OO 

i f , A inh A Vt+ {0) exp{i(ip - vtx)} C exp{z^(f - f p )} 

2^2 J av [^ + + V + i leg y + AU J [a + {u) + a_(fi) - 2tSk] 

— oo — OO 

Q °ri 7e9 t A \"^ {exp[-h<*+(")z + *Skz] - exp[-i a+ (u)L + ±a(Q){z - L) + tSkL}}, (59) 

where u;_ = ui eg + rj(uj eg — uj + ) is the carrier frequency of the echo held, ip = (l + r))(uj e g — u>^-)ti, fli (zv) = (u;_ — uj eg ) + 
i](lu + — uj eg + v) + i(l + ri)"f eg , and a_ (u) = r- + ^A°^I+i(^'"-L +v)] IS * ne absorption coefficient of the atomic system 
during echo emission. It is characterized by the compressed inhomogeneous broadening ryAj„/j. Assuming again the 
inhomogeneous broadening to be large compared with the initial light spectrum, Ai n h 3> Sujf, we can simplify the 
dependence on the phase mismatch factor: 

2 ^2 [7 eg + xAinh + i(u eg + 0)] [^ eg + A lnh + i(uj eg - u + - v)] 



a+{v) + a_(fi) - 2iSk a a { leg + A inh ) [ 2 j eg + (1 + X )A lnh + i(2uj eg - cj+ - cj_ - z/ + ft) ~ 2i8k {rleg+ * Ainh) ] ' 

(60) 



The main part of the integration in Eq. ( 59 1 over in the complex plane is determined by the pole singularity 
at — Qi(u), while the other singularities describe fast decaying signals associated with the large inhomogeneous 
broadening. After the integration we obtain the following echo field amplitude at the medium output [z = 0): 



A-(f,z = 0) = 



r]exp{iip - (1 + vhegjr - h)} 
7T [(77 + 1) + (1 - ri) leg /A inh ] 



exp{w[7?(f -h) -exp[-|a(i/, rf)L + iSkL}} - 

(r ) +l) + (l-r 7 ) 7e9 /A I „ fe ZZ (r ) +l) + (l-r ) ) 7e9 /A, rlh / 

where a{y,rj) — [a + {i>) + a_(f)] describes the combined influence of absorption and dispersion effects on the echo 
field. 



Eq. (61) allows calculating the compression efficiency. For instance, setting 77 = 1, we can derive the echo field 
in the presence of phase mismatch [35l [36] , which limits the recall efficiency. Setting 8k = 0, we find the echo to be 
a perfect, temporally reversed copy of the input field, provided the optical depth is sufficiently large for all spectral 
component of the field (a.(y,rj)L ^> 1. This confirms previously obtained results [HI 1551 157] . For sufficiently large 



inhomogeneous broadening Sujf/A in h >lwe find a(v, rf) = (1 + X/rf)a and Eq. (61 1 can be simplified to: 



A-{f,z = 0) = 
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Sofo^M a a L) exp{iip - (1 + ?7)7e S (f - t 1 )}^A+[-» 7 (f - h) + tj, 



(62) 



where M^> (8k, a a L) and e\/ 2 (ri) are given in Eqs. (32 ) and (33 ). While we are only interested in the quantum state of 



the retrieved and compressed photonic qubit (as described by Eq. (28)), which we obtain through Eq. (62 1 assuming 



initial photonic time-bin qubits (Eq. \\2\ ) with Gaussian shape of the basis wave packets, we note that one can also 
derive the amplitude of the atomic state (Eq. (24)) after the echo signal emission by inserting Eq. (59) into Eq. (26). 



B. Appendix - Longitudinal broadening 



As in the above discussed case of transverse broadening, we are only interested in the retrieved photonic time-bin 
qubit state. Averaging Eq. (26) over the phase fluctuations, and after a variable transformation to a moving reference 
frame r = t— z/c, z — z, the light-atom equations for t < t\ become 



J^i(t, z) = i(irn Sg*/c)b (T, z), 



(63) 



:b (r,z) = (iyz -% g )b (T,z) +igAi(r,z), 



(64) 



where we used the substitution bj(r) — b (r, Zj) exp{— iu)2iT} with b (r, z) describing averaged local atomic coherence, 
and Ai(t, z) = A + (t, z) exp{— i(uj+ — w eff )r}. The field amplitude A + (t, z) is given by Eq. ( 12 ) with additional factor 
exp{i4>i g } where 4>lg = \u+L/c. 



The general solution of Eqs. (63) and ( |64[ ) is obtained using a temporal Laplace transformation jJU]. We are 
interested only in the atomic coherence excited by the input light, which gives rise to the echo field irradiated after 
further evolution, while the initial light field disappears in the medium (A±(t ^> St, z) — \ 0): 



b (n,z) = - 



2?r 



duje 



{0J + X Z + ijeg) 1+lC/x 



Ax{u,-L/2). 



(65) 



Ai(uj,~- L/2) = f x ' oo dTe lu ' T Ai(T:—L/2) is the field spectrum at the input of the atomic medium (z — —L/2), and 
C = ■nn Sg 2 / c. Similar to the transverse scheme, we find the atomic coherence on the resonant optical transition 1 — > 2 
at the moment of re-emission of light using the phase mismatch factor Skz: b (ri + (e — > 0)) = 6 (ri) exp{iSkzj}. 

We now analyze the backward echo emission. After changing the atomic detunings A(t > t\) — x'z, we find the 
following system of equations for the atomic coherence and irradiated echo field A 2 (t > t 1; z) (where f = t + z/c): 



-§ i A 2 (f,z) = -i(irn Sg* /c)b (f,z), 



(66) 



■§fb {f. z) = -{ixz + J eg )b (T, z) + igA 2 (f, z) 



(67) 



where A 2 (t, z) = A_(t, z) exp{— i(uj_ — uj eg )f}, and w_ is the new carrier frequency. Using a Laplace transformation 
in the solutions of Eqs. (66) and ([67]), we find Aiij , z < L/2) as a function of the stored atomic coherence 6 (fi, z): 



A 2 {t,z = -L/2) 



C 



L/2 



{2-Kf 



dz 



-L/2 



i { -x'L/2 - m - i leg ) V 

(X'Z 1 - U) - ijeg)^'^ 



exjp{iSkz'} 



dio'e 



otl (^- X L/2 + z 7eg )'C/x 



(68) 
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where we used n = tx. Taking into account backward emission of the echo field, as well as large inhomogeneous 
broadening with respect to the input light spectral width, and using the simplifications 

(a/ - xL/2 + i leg ) K/x = (xi/2) lC/x exp{ 7es r m - ir(/ X } exp{-n- m u/}, 
(-x'L/2 - w - i7e9 )- 4C/x ' = MLM-KM exp{ 7eg r^ - exp{-fr», 



where r' m = (C/xO [( X 'L/2)2+ 7 g ] - and r ™ = (C/x) [( X L/2)^+ 7 g ] - fe^J' we find > aftcr chan S in g the order 01 



(C/x') 



(C/x) 



the integrations in Eq. ( 68 ) from tj, z', ui' to u> f , z' , ui: 



A 2 (r,z = -L/2) = ^ expfr^ + r m ) <(!/*' + 1/ X )}(x' L/2)-^' ( x L/2)^ 



du/exp{-m/(ii + T m )}Ax(u}', -L/2) 



L/2 



dz - 



cxp{iSkz'} 



dus 



exp{-io;(r -tx + T' m )} 



-L/2 



{(J + X Z' + ileg) 1 + K/x J {X'Z' ~ W ~ *7e») 1_<C/x ' 



Taking again into account that X L 3> Sujf, and using the tabled integral 



(69) 



{^L/2)-' K/x ' / du 



exp{mT} 



2n Vx (T) 



- {T X 'L/2)- lC/x ' exp{- 7e9 T + \^/x'}, 



(C/x')n-K/x 



where r\ x (T > 0) = 1 and rj x (T < 0) = is a Heaviside function, we integrate Eq.(69) over ui, then over z' and a/, 
leading to 



A 2 (t,z = -L/2) 



2^x , exp{-^7^C/x , - X/x> 

C r[-zc/x']r[zc/x] 



t] x (t -tx + T m - Sk/x') 



(( X 'L/2)(r - tl + r' m ))'^' (( X 'L/2)\r - tx + r' m - 6k/ X '\) lC/x 



exp{- 7es (l + t?)(t - ti) + leg r]T z }Ax{-r)[r - (1 + - r 2 ]}, (70) 

where r\ = x'/x is the compression factor, and F[±ix] are Gamma functions. Interestingly, as described by Eq. ( |70[ ), 
phase mismatch does not effect to the quantum efficiency of the echo emission but only leads to a temporal shift of 
the echo emission to T' echo + t z (where T' echo = (1 + l/rj)tx, t z — Sk/x' + T m/v — T' m ). However, this shift results in 
an additional nonlinear phase shift due to the deviation from perfect temporal reversibility. Furthermore, we find 
that the echo emission is conditioned on a phase mismatch Sk/x' > T L + $t — (*i + r m)/?7- Taking into account that 
Tm < St and r' m < St, we find Sk/ X > —tx- 

Finally taking into account the relations between ^(r,^) and A-(f,z), and between Ax{f,z) and A + (f,z), we 



find the new carrier frequency cj_ 



r?(o 



oj + ) of the echo signal, which coincides with the case of transverse 



broadening. We also find A_(f, z) — A2(f, z) exp{iri(u) eg — lo + )t}. Using this relation, together with Ai{—t](t — 
T echo - T z)} = A+{-ti(t - r' echo - t z )} exp{z?7(w + - w eff )(r - r ' h - t z )} in Eq. ^70| ), we find the solution for A-(f, z) 
as a function of A + . Putting the obtained solution into Eq.(ll), assuming again a Gaussian shape of the photonic 



wave packets, we obtain the solution for the retrieved photonic time-bin qubit state given in Eq. (40 1 after a simple 
algebraic calculation. 



